Stability of Large Flocks: an Example 



J. J. P. Veerman*(l,2) 
(1) Dept. of Math. & Stat., Portland State University, Portland, OR 97201, USA. 
(2) Center for Physics and Biology, Rockefeller University, NY, NY 10021. 

February 9, 2010 



Abstract 

The movement of a flock with a single leader (and a directed path from it to every agent) can be stabilized. 
Nonetheless for large flocks perturbations in the movement of the leader may grow to a considerable size as they 
propagate throughout the flock and before they die out over time. As an example we consider a string of TV + f 
oscillators moving in IR. Each one 'observes' the relative velocity and position of only its nearest neighbors. This 
information is then used to determine its own acceleration. Now we fix all parameters except the number of 
oscillators. We then show (within a certain class of systems) that a perturbation in the leader's orbit is almost 
always amplified exponentially in N as it propagates towards the outlying members of the flock. The only exception 
is when there is a symmetry present in the interaction: in that case the growth of the perturbation is linear in N. 

1 Introduction 

We take up the study of a system very familiar from physics ([5] for example): a long string of coupled linear 
oscillators (called 'agents' or 'cars' from here on) with linear nearest neighbor interactions. In order to model things 
like traffic we change certain features of the model that is so familiar in the context of statistical physics. The 
most significant ones are that we do not assume that interaction is symmetric and that the oscillators are damped. 
Furthermore we allow the system to have non-zero mean velocity. Finally we assume that the chain is finite and we 
take into account boundary effects from the ends of the chain (no periodic boundary conditions). The class of systems 
whose study we take up in this paper is illustrated in Figure 11.11 and made precise in Equation 12.11 
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Figure 1.1: In the upper figure the 'symmetric' (p — \) flock is illustrated. Each agent is linearly coupled to its nearest 
neighbor. At t = the agent labeled 0, the leader, undergoes a forced motion: either an oscillation or a kick in the 
direction of the arrow above it. It receives no feedback from the flock. The asymmetric interaction is suggested in the 
lower figure, where we drew the k-agent with its interactions and their weights. The arrows give the direction of the 
information flow. 

The purpose here is to gain a qualitative understanding of how a perturbation in the movement of one the agents 
(which we call the leader) propagates throughout the flock of agents. Here is a more precise formulation. Suppose 
that a flock is moving coherently or "in formation" , ie: each member has the same constant velocity, and inter-agent 
distances are constant. Now the leader's motion is perturbed. After a while also the agent farthest away from the 
leader (called the trailing car) senses the effect, and its orbit will start to deviate from the in formation orbit. What is 
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the ratio between these two perturbations as a function of the size N of the flock, while keeping all other parameters 
fixed ? 

The immediate motivation is to study which kinds of interactions might take place in actual large flocks or 
might seem desirable to implement in artificial flocks (think of cars on a highway equipped with an automatic pilot). 
It would seem that in both types of systems the kinds of interaction preferred would be those where the growth rate 
just mentioned is as low as possible. 

Ideally we find the growth rate of perturbations in many different models, for example where interactions are 
allowed to include up to 2 or 3 or 10 neighbors on either side. This seems analytically out of reach. Thus the 
acceleration of the k-th car is determined by its observations of the k — 1-st and the k + 1-st car. What we chose to 
do instead is to introduce a weighted nearest neighbor interaction, that is: For p £ [0, 1] we weight the information 
coming from the k — 1-st car with 1 — p, and the other with p (see Equation (|2.1[1 ). 

This line of thought has many potential applications in technology. Suppose for example one wants to automat- 
ically control dense traffic on a single lane road. Each car is equipped with sensors that register the relative velocity 
and position of its neighbors. This information is then used to regulate the acceleration of the cars. A typical example 
is the 'canonical traffic problem' proposed in [7] and [8], where the lead car abruptly accelerates (when, for example, 
a traffic light turns green), and the other cars try to follow. 

Our conclusions for the systems typified by Equation l2.ll can be summarized as follows. In all systems the ratio 
of the size of perturbation of the trailing car to the perturbation of the leader grows exponentially in the size TV of 
the flock, with one exception: when p—\- In the last case (called the symmetric case because equal attention is paid 
to the cars on either side) this growth is only linear as function of N. The study of the symmetric case was reported 
in [318]. Here we continue that program by investigating the asymmetric case. 

The curious fact that perturbations in the leader's position or velocity is necessarily unbounded as the size of 
the flock grows is of course of paramount importance in many applications (automated traffic, biological flocks, and 
so on) . In fact this has been commented upon by several authors already, notably [H [5] in cases similar to our p = 
and p = \ cases. In this note we give precise estimates for how these perturbations grow as function of the number of 
agents, not only in those cases but also for other models (all p £ [0, 1]). 

A more distant motivation to study these systems is to attain a general understanding of linear oscillators 
interacting according to a more general graph (called the "communication graph" , see for instance [B] for definitions) . 
In that context similar questions arise but in a more general context. What we present here is a simple example. 

Section [2] defines the model. In Section [3] the definitions of stability we use are given and briefly discussed. In it 
we also state the main result of this note. Sections [4] and [5] prove these results. (Some of the more calculational steps 
in these proofs are relegated to the Appendix.) 

Notational Conventions: To avoid confusion, we finally list two important conventions here. The first is that after 
Section 3] we assume that both / and g are negative reals to insure asymptotic stability (Theorem I4.2[) . The second 
is that in order for certain expressions {p+{iuj) and fj,— (iu>), for uj > 0) in Theorem 15. 2p to be continuous functions it 
is convenient to define the symbol *Jz as the root with angle in the interval [0, it) (branch cut along the positive real 
axis). 
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2 The Model 

We define the model and propose a notion of stability for flocks. Our strategy is to study qualitative aspects of the 
solution for fixed p, /, and g, and as we let N tend to infinity. 
Let / and g be real, and p £ [0, 1]. The model is given by: 

Vi £ {1, • • • N — 1} : Xi = m 

Ui = f {{xi - hi) - (1 - p){xi-\ - hi-x) - p{x l+1 - h i+ x)y +g{ui - (1 - p)ui-x - pu i+ {\ 
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and 



x N = u N 

un = f {{xn - h N ) - (xjv-i - hn-x)} + g{u N - ujv-i} 
x Q = x (t) given 



(2.1) 



The (constant) parameters hi determine the desired relative distances between agents i and i — 1 as hi — hi-i. The 
feedback parameters / and g are independent of i and time. Note that the total number of agents is in fact also a 
parameter in this problem. We do not carry this into the notation. The subscript N will always stand for the last 
agent in a system with agents numbered from to N. Similarly we do not carry the dependence on the parameters / 
and g explicitly into the notation. 

One can show that the orbits of this system such that the distances between successive cars are preserved 
(namely: — (hi — /ij-i)) form a 2-parameter family, namely Xk (t) = xq(0) + Vo(0)t + hk andifc(t) = i>o(0). We will call 
these orbits in formation orbits. 

It is advantageous to write Equation (|2.1j) in a more compact form. Introduce the notation 



z = (zi, ii, z 2 , z 2 , ■ ■ ■ , z N , z N ) = (x% - hi, ui, x 2 - h 2 , u 2 , ■ ■ ■ , x N - h N , u N ) 

The leading car is not encoded since its orbit is a priori given. The system can now be recast as a first order ODE: 

z = Mz + T {t) . (2.2) 

The details of this are discussed in [6], here we just give the relevant definitions. 

Let / and P arc iV-dimcnsional square matrices, where / is the identity and P is given by 



P 



where Q t - 
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V 



\ 



p 
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(2.3) 



P = I — Q p is called the reduced graph Laplacian. It describes the flow of information among the agents, with the 
exception of the leader (hence the word 'reduced'). The 2x2 matrices A and K are given by: 



1 




and K 





/ g 



(2.4) 



The orbit of the leader is assumed to be a priori given and therefore only appears in the forcing term To(i). We will 
refer to this agent as an (independent) leader. Analyzing Equation (|2.1[) and assuming without loss of generality that 
ho = 0. one gathers that: 

/ \ 

(1-p) (fz (t)+gz (t)) 
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To define M of Equation (12.21) in terms of these quantities, we use the Kronecker product ((g)) 

M = I ® A + P ® K . 



(2.5) 



(2.6) 



The advantage of this somewhat roundabout way of defining the matrix M is that in the eigenvalues of the reduced 
Laplacian P in many cases are known. From that the eigenvalues of M can then be derived. That is the program 
followed in the Section @] 



3 Stability of Flocks 

Definition 3.1 The system given in Equation \2.2\ is called 'asymptotically stable' if all eigenvalues of M have negative 
real part. 
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Suppose for the moment that To(i) — fot t > to. Then the solution of the system tends to exponentially fast (in 
t) if and only if the system is asymptotically stable. This corresponds to the usual notion of asymptotic stability (see 
for example [I], Section 23). The question when the system is asymptotically stable has a straightforward answer (see 
Theorem 14. 2[) : it is if and only if / and g in Equation (|2.ip are negative. We will therefore from now on assume that 
f and g are negative. 

One can show (Proposition I5.1|) that if the leader executes an oscillation of the form e lut , then Zk(t) tends to 
afe(iw)e laj * as t tends to infinity. The functions ak(iu>) are called the frequency response functions. 

Definition 3.2 Let An = sup wg jj> |ajv(*w)|- The system is called 'harmonically stable' if it is asymptotically stable 
and if limsupjy_>.oo I^ljvl 1 ^ < 1- Otherwise the system is called 'harmonically unstable'. 

This kind of instability roughly says that certain long-term oscillatory perturbations in the orbit of the leader will 
have their amplitude magnified by a factor that is exponentially large in N. In Theorem 15.51 wc establish that the 
system is harmonically unstable if p ^ 1/2; harmonic stability for p = | was established in [7]. 

In earlier work ([5]) we studied a 'fundamental traffic problem' which roughly corresponds to setting the accel- 
eration of the leader equal to the Dirac delta function, 6(t), and zq(0) = zq(Q) = 0. This gives Zo(t) = t for t > 
which can be substituted into Equation (|2.5|) . 



Definition 3.3 Consider Equation 12. 2\) with forcing determined by zo(t) = S(t) and subject to the initial conditions 
Zfc(0) = ifc(0) = and ifc(0) = 0. Let Z^ = sup t>0 \~(z^(t) — Zo(t))\. The system is called 'impulse stable' if it is 



asymptotically stable and if for i 0, 1 and 2, we have \\msn~PN^ roo Z^ < 1. (And 'impulse unstable' in the other 
case.) 

Loosely interpreted this kind of instability means that if we give the leader a 'unit- kick', then that perturbation travels 
through the flock and causes sup |xjv(t)|, sup |iAr(i)|, or sup |xjv(i)l *° grow exponentially in N, before eventually dying 
out (due to asymptotic stability). 

Impulse stability is perhaps at first sight more natural or appealing than harmonic stability because it is formu- 
lated in the time domain, whereas the latter takes place in the Fourier domain. However mathematically the criterion 
is much harder to check. In [8] we proved that the case p = 1/2 is impulse stable. But the case p / 1/2 is much more 
problematic and will be taken up in a separate work ( [5] ) . 

It therefore may be argued that all three kinds of stability are necessary to form large flocks. From the above 
remarks, it follows that of the systems investigated here only the one with p = \ satisfies all three. It is interesting 

that even in that case linear growth of zff still takes place (see [TJIH])- This seems to be the best case possible, we 
summarize this discussion with our main result. 

Theorem 3.4 The system given by Equation h2.1\) satisfies all three stability criteria if and only if f and g are 
negative and p = 1/2. 



l/N 



4 Asymptotic Stability 

We derive the criterion for asymptotical stability. 

In the statement of the next result we use the following equation, where p G (0, 1) and <j> are real variables: 

(2p- l)cot<£ = cot Ncj) . (4.1) 

Recall that the matrix P is defined in Equation (|2.3[) . 



Proposition 4.1 (J^j) For any p <G (0, 1), the matrix P has N distinct eigenvalues {A^j^Q 1 : 

i) If pe (0,i]= fori£{0,...,N-l}, X t = l-2y/p(l-p) cos <f> t , where fa e (f ,^ L ) solves (£7J). 

ii) If p G (h, ^ff-]'- Identical to i). 
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iii) If p G (^±1,1): for I G {1,...,N- 2}, \ £ = 1 - 2y/p(l-p) cos0 £; wfcere 0, G (f , sotoes |£7P; 

Ao = ^ (V) ^ + (( 1 7 £ ) 2Ar_2 ) aW A ^ = 2 - Ao - 

One can show (see [BJ [71 [5]) that the eigenvalues of M defined in Equation (|2.2[) are given by the solutions Vt± of 

i/ t -X t gv-X t f = , (4.2) 

where A^ runs through the spectrum of P. So: 
Theorem 4.2 The eigenvalues of M are 




where Xp runs through the spectrum of P. Because the Xi are contained in the interval [0,2] (see Proposition ^. 1\ ) , the 
system is stabilized (or asymptotically stable) if and only if both f and g are strictly smaller than zero. 



5 Harmonic Stability 



In this section we calculate (Proposition 15. 2[) and study the properties (Theorem I5.5[) of the frequency response 
function of the trailing car in detail. The main result of this Section is Theorem 15.51 that says that if [0, 1)\{±} then 
the system is harmonically unstable. (If p = 1 the question is moot as the leader's motion goes unperceived by the 
flock.) 

The following constant will frequently simplify formulae: 

_l-p 1 

k = or p — 

p 1 + K 



Proposition 5.1 If zo(t) = e vt and v € iM and v spec (M), there is a sequence a — (ai ■ • ■ , ajv) of complex 
numbers {ctk( l ')}k = i so that trajectory of the system is asymptotic to: 

Zk(t) — a k (v) e vt where a(i>) = — (M — vl)~ 1 g . 



Proof: (See also [BJ.) Since Zo(t) = e vt , the non-autonomous term in Equation (|2.2[) can be written as To(t) = goe vt , 
where go is a constant vector. The general solution of the system is: e Mt zo + (M — vl)^ 1 (e Mt — e vt \ go, and e Mt 
tends to zero. ■ 



The complex functions ak{v) (where v G i]R) are called the frequency response (of the fc-th agent). 

Proposition 5.2 i): For p £ (0, 1)\{^} the frequency response function of the k-th agent is given by the functions: 



M+ - M+ ) M+ - (M- - M- ) M- 

where p± = p±(v) = ^- (7 ± ^ j 2 - 4p(l - pf) and 7 = 7(1/) = - + 9U — — 

2/0 V / / + gv 
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(W$) When p = h the above expressions simplify to: 



a k\ v ) = — 7WZ — N— 



f + gv- v l 



where P±M — 7 ± vT-l and 7 = jM 

f + gv 

(J7j) When p = 0: 



Proof: The reasoning of i) is identical to that in Lemma 3.2 of [5]. ■ 

Remark: Even though p± are not rational functions of v, one can check that in fact the a k {v) in fact are proper 
rational. 

Remark: The coefficients p± in the Theorem are the roots of the following quadratic equation: 

pp 2 - 1P +(l- P ) = Q . 



The most important case of Theorem 15.21 is the frequency response of the trailing car, labeled N , when p € 
(0, 1)\{^}. The previous result immediately implies: 

Corollary 5.3 For p <E (0, the frequency response function of the last agent is given by 

1 + K N p + — H- 

a*N\ v ) = K 



(m+ - M+ - (M- - M- 1 ) 
with p± and 7 as before. 

Because the inverse Fourier Transform of a at equals the real valued function z/v(i), we are mainly interested in 
the case where v = iuj and w is a positive real number (the frequency). Note that then also q.at(— iu>) must equal the 
complex conjugate of a^(ioj). In fact, using Lemma |6. II in the previous Proposition yields that ajv(0) = 1. Thus it is 
sufficient to study a^{iui) only for to > 0. 

Proposition 5.4 Suppose f , g, lu > and p £ (0, 1/2) U (1/2, 1) are all fixed. Then, for r £ (0, 1) as in Lemma Wlk 
as N large tends to infinity: 

a N (iw) = p_ ^(1 + Olr )) . 

Proof: Use Proposition 15.21 and the fact that p~p+ = K to rewrite 

M^) = — 1 zr — ■ (5.1) 

Since 

/it- - /il 1 -1 - ^ 2 

it suffices to prove that if w > 0, then p + (uj) 2 7^ 1. Now suppose that p + (uj) 2 = 1, then the second remark after 
Proposition 15.21 implies that 7(iw) = ±1 and so Lemma RT21 implies that w = 0. ■ 

Remark: It is important to realize that the even for large but fixed N , the last factor could become large if we let 
/U_l_ approximate 1. The results in the appendix can be used to show that this can happen if and only if p < j and uj 
very close to zero. We make use of this fact in the proof (for p < i) of the next and main result of this section. 
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Theorem 5.5 For all p £ [0, 1)\{|}, ^jv grows exponentially in N. When p = i, growth is linear in N. 

Proof: The second statement has been proved in [8]. Fix p £ (0, 1/2) and let lu + as in Equation 16. II Lemma 16741 and 
the remark thereafter now imply that \p-(iu>)\ > 1 if and only if ui £ (0, uj+). The result follows directly from the first 
part of Proposition 15 .41 Finally fix p £ (-|, 1), or k £ (0, 1). First use Lemma to see that if w 2 = ^|/|(1 — k) 2 k" _1 , 

then n + = 1- {1 ~f ] k n - x +0(k 3N / 2 ) and p_ = k(1+ {1 ^ 2) k n - 1 ) + 0(k 3N / 2 ). Substitute this into the denominator 
of on in Corollary 15.31 The leading order cancels. The next term is of order at least k 3N / 2 . Thus An is of order at 
least k~ n / 2 . ■ 



6 Appendix: Technical Results 



For completeness we collect a number of straightforward results that are necessary for development of the theory, 
but would clutter the exposition in the main text. Various relevant quantities are evaluated for v — iuj where uj is real 
and non- negative. As observed in the main text, we use to > without loss of generality. The conventions mentioned 
in the Introduction also hold. 



Lemma 6.1 i): For p £ (0, |) and lu > small: 
ii): For p £ (i, 1) and ui > small: 



Proof: By sheer calculation. (See [7] for some of the computational details.) 
Remark: This expansion diverges for p = 1/2; in that case we have ([7]): 



Lemma 6.2 7 iw = 1 - — + z — 

/- + w V /- + w V 



Lemma 6.3 For eac/i p G (0, 1)\{|}, the number r = sup w exists and is contained in (0, 1). (See Fiaure [6J\ ) 

Proof: From Lemma HT2l j(iuj) ~ — -y when w is large. Substitute this into the expression for p± in Theorcm l5.2l to 

see that for large ui, in fact r^T^ becomes very small. When u = 0, Lemma l6Tl implies that r^l^r! = minjK, k -1 j. 

It is now sufficient to prove that for uj £ M + the absolute values \p± \ are never equal. So suppose there are ojq 
and 9 £ IR so that p+(iuo) — p-(iojo)e ie = 0. The first item in Theorem 15.21 gives : 

7(l-^) = -V / 7 2 -4p(l-p)(l + e 4e ) . 
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Figure 6.1: The eigenvalues p, + {iuj) (blue) and H-{iuj) (red) of C when f = g = — 1 for uj positive. From left to right: 
p = 0.4, 0.5, and 0.6. In addition the circles with radii and n are drawn in green and black, resp., where k = — = £. 



Dividing this by 1 + e l6 ', squaring the equation, and noting that J^. = — (tan(|)) 2 , we see that 



I tan 5 



4p(l - p) 



This implies that 7 is a positive real and therefore 7 is real for some uj 7^ 0, which is impossible by Lemma 16.2 



Lemma 6.4 For each p £ (0, 1/2), there is a unigue uj + > such that 

uj € (0,uj + ) =£> > 1 and uj > uj + =>■ \p,^(iuj)\ < 1 



Proof: We know that p.-(0) = 1 and (from the proof of the previous Lemma) for large uj: \/j,-(uj)\ is small. It is 
sufficient to prove that uj+ is the unique solution in (0, 00) of \p,-(iuj)\ = 1 and that it is simple. 

Consider the characteristic equation pp 2 — 7/i + (1 — p) = and suppose that there is a root p = e lfl '. Then 
7 = pe 10 + (1 — p)e~ %e = cos(#) + i(2p — 1) sin(#). Equate this to the expression given in Lemma IB~2l and use the fact 
that cos 2 (#) + sin 2 (6>) = 1 to obtain: 



1 - 



^ 2 |/l 



f 2 +uj 2 g 2 J (2p-l) 2 \f 2 +uj 2 g 2 



= 1 



This equation factors as follows: 



UJ' 



■^ 4 + (/ 2 -2|/| 5 V-2|jf =0 



,(2p-l) 2 

The second factor gives exactly one simple positive root for uj 2 , yielding a unique simple positive root uj = uj + . 



Remark: In fact. 



uj\ = (l-2p)|/| 1 



2g 2 



(1 - 2p) + J 1 



l/l 



9 2 



(6.1) 
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